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Abstract.  In  this  paper,  Bezier  curves  and  surfaces  have  been  constructed  based  on 
modified  Bernstein  bases  functions  with  shifted  knots  for  t  €  [yqrg,  •  Various  prop¬ 
erties  of  these  modified  Bernstein  bases  are  studied.  A  de  Casteljau  type  algorithm  has 
been  developed  to  compute  Bezier  curves  and  surfaces  with  shifted  knots.  Furthermore, 
some  fundamental  properties  of  Bezier  curves  and  surfaces  with  modified  Bernstein  bases 
are  also  discussed.  Introduction  of  parameters  a  and  /3  enable  us  to  shift  Bernstein  bases 
functions  over  subintervals  of  [0,1].  These  new  curves  have  some  properties  similar  to 
classical  Bezier  curves.  We  get  Bezier  curves  defined  on  [0, 1]  when  we  set  the  parameters 
a,/3  to  the  value  0.  Simulation  study  is  performed  through  MATLAB  R2010a.  It  has 
been  concluded  that  Bezier  curves  that  are  generated  over  any  subinterval  of  [0, 1]  based 
on  modified  Bernstein  bases  functions  are  similar  to  the  Bezier  curves  that  are  generated 
based  on  classical  Bernstein  bases  functions  over  the  interval  [0, 1]. 
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1.  Introduction  and  preliminaries 

It  was  S.N.  Bernstein  [1]  in  1912,  who  first  introduced  his  famous  operators 
Bn  :  C[0, 1]  — >  C[0, 1]  defined  for  any  n  E  N  and  for  any  function  /  e  C[ 0, 1]  where 
C[0, 1]  denote  the  set  of  all  continuous  functions  on  [0, 1]  which  is  equipped  with 
sup-norm  ||.||C[o,i] 

B»(/;x)  =  y  (  x  6  [0, 1],  (1) 

k= 0  '  '  '  ' 
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and  named  them  Bernstein  polynomials  to  prove  the  Weierstrass  theorem  [9]. 
Bernstein  showed  that  if  /  £  C[0,  1],  then  Bn(f-,x )  ^4  /(x)  where  ”  ^4  ”  repre¬ 
sents  the  uniform  convergence.  One  can  find  a  detailed  information  about  the 
Bernstein  polynomials  in  [10]. 

Later  it  was  found  that  Bernstein  polynomials  possess  many  remarkable  prop¬ 
erties  and  have  various  applications  in  areas  such  as  approximation  theory  [9], 
numerical  analysis,  computer-aided  geometric  design,  and  theory  of  differential 
equations  due  to  their  fine  properties  of  approximation  [23]. 

In  computer-aided  geometric  design  (CAGD),  Bernstein  polynomials  and  their 
variants  are  used  in  order  to  preserve  the  shape  of  the  curves  or  surfaces.  One 
of  the  most  important  curve  in  CAGD  [27]  is  the  classical  Bezier  curve  [2]  con¬ 
structed  with  the  help  of  Bernstein  basis  functions.  Other  works  related  to  differ¬ 
ent  generalizations  of  Bernstein  polynomials  and  Bezier  curves  and  surfaces  can 
be  found  in  [3,  4,  5,  7,  8,  11,  12,  13,  18,  20,  21,  23,  24,  26] 

A  new  analogue  of  Bernstein  operators  using  the  concept  of  post  quantum 
calculus  ((p,  g)-calculus)  in  approximation  theory  have  been  introduced  recently 
by  Mursaleen  et  al.  in  [12].  Later,  based  on  (p,q)- integers,  some  approximation 
results  for  Bernstein-Stancu  operators,  Bernstein-Kantorovich  operators,  ( p ,  q )- 
Lorentz  operators,  Bleimann-Butzer  and  Hahn  operators  and  Bernstein-Shurer 
operators  etc.  have  also  been  introduced  in  [13,  14,  15,  16,  17]. 

Also  see  a  recent  work  in  approximation  theory  On  Approximation  by  Stancu 
type  Jakimovski-Leviatan-Durrmeyer  operators  [19]. 

In  [7,  8],  Khalid  et  al.  have  shown  applications  of  post  quantum  calculus  in 
computer-aided  geometric  design  in  terms  of  flexibility  and  applied  these  Bern¬ 
stein  bases  for  construction  of  ( p ,  q)- Bezier  curves  and  surfaces  which  is  further 
generalization  of  g-Bezier  curves  and  surfaces  [21,  23].  For  other  relevant  works 
based  on  Bezier  curves,  we  refer  the  readers  to  [4,  5,  11,  18,  21,  22,  23,  25]. 

In  1968,  Stancu  [28]  showed  that  the  polynomials 


(d“'®/)M  =  ±  (")** (1  -  x)-‘/(D t|),  (2) 

converge  to  continuous  function  /(x)  uniformly  in  [0,1]  for  each  real  ck,  / 3  such 
that  0  <  a  <  (3.  The  polynomials  (2)  are  called  Bernstein-Stancu  polynomials. 

In  2010,  Gadjiev  and  Ghorbanalizadeh  [6]  introduced  the  following  construe- 
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tion  of  Bernstein-Stancu  type  polynomials  with  shifted  knots: 


Sn,a,p{f;x)  = 


n  +  fi2 


n 


n  n 

E 

k= 0 


x  — 


0-2 


n  +  a2 


n  +  /32J  \n  +  (32 


—  X 


n—k 


k  +  ai 

n  +  j3\ 


(3) 

where  nffip2  <x<  r^a22  and  ak,fik  (k  =  1,2)  are  positive  real  numbers  provided 
0  <  a\  <  a2  <  Pi  <  (32.  It  is  clear  that  for  a2  =  fi2  =  0,  the  polynomials  (3) 
turn  into  the  Bernstein-Stancu  polynomials  (2)  and  if  a\  =  a2  =  /3i  =  f32  =  0, 
then  these  polynomials  turn  into  the  classical  Bernstein  polynomials. 

In  recent  years,  generalization  of  the  Bezier  curve  with  shape  parameters  has 
received  continuous  attention.  Several  authors  were  concerned  with  the  problem 
of  changing  the  shape  of  curves  and  surfaces,  while  keeping  the  control  polygon 
unchanged  and  thus  they  generalized  the  Bezier  curves  in  [7,  8,  22,  23]. 

The  outline  of  this  paper  are  as  follows.  Section  2  introduces  a  modified 
Bernstein  functions  with  shifted  knots  a  ^  and  their  properties.  In  Section 
3,  Bezier  curves  based  on  modified  Bernstein  bases  alongwith  degree  elevation 
and  a  de  Casteljau  algorithm  are  presented.  In  Section  4,  we  define  a  tensor 
product  patch  based  on  Algorithm  1  and  we  discuss  its  geometric  properties 
as  well  as  a  degree  elevation  technique.  Furthermore,  tensor  product  of  Bezier 
surfaces  on  ^§]  x  [^qrg,  f°r  Bernstein  polynomials  with  shifted  knots 

is  introduced  and  its  properties  inherited  from  the  univariate  case  are  discussed. 
In  Sections  5  and  6,  we  have  given  concluding  remarks  and  MATLAB  codes. 

In  next  section,  we  construct  basis  functions  with  shifted  knots  with  the  help 
of  (3). 


2.  Modified  Bernstein  bases  functions  with  shifted  knots 


The  modified  Bernstein  bases  functions  with  shifted  knots  are  defined  with 
the  help  of  (3)  as  follows: 


G 


k 

n,a,/3 


(  a  \A/n  +  a  \ 
\  n  +  /3 )  \n  +  /3  / 


(4) 


where  ^  <  t  <  and  a,/3 


are  positive  real  numbers  provided  0  <  a  <  f3. 


2.1.  Properties 

Theorem  1.  The  modified  Bernstein  bases  functions  with  shifted  knots  possess 
the  following  properties: 

1.  Non-negativity:  Gkn  a^  (t)>  0  k  =  0, 1,  •  •  •  ,  n,  t  e  [n'f3. 
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2.  Partition  of  unity : 


(<)  =  !.  for  every  t 


k= 0 


3.  End-point  interpolation: 


Gk 
n,ot,/3 


a 


n  +  f3 


a  n  +  a 


n  +  f3  n  +  (3 


1,  if  k  =  0 
0,  k  +  0 


Gt 


n  +  a 
n  +  /3 


1, 


if  k  =  n 
k  /  n. 


4.  Reducibility :  when  a  =  (3  =  0,  formula  (4)  reduces  to  the  classical  Bernstein 
basis  on  [0, 1]. 


Proof.  All  these  properties  can  be  deduced  from  (4).  ◄ 


Figure  1:  Modified  Bernstein  bases  functions  with  shifted  knots 

Figure  1  shows  the  modified  Bernstein  bases  functions  of  degree  3  with  shifted 
knots  for  a  =  4,  (5  =  6.  Here  we  can  observe  that  the  sum  of  blending  functions  is 
always  unity  and  also  satisfies  end  point  interpolation  property.  In  case  a  =  /3  = 
0,  it  turns  out  to  be  classical  Bernstein  basis  on  [0, 1]  which  is  shown  in  Figure  2. 
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Figure  2:  Classical  cubic  Bernstein  bases  functions 


2.2.  Identities 

Below  we  state  some  important  identities  which  we  will  use  later. 


Proof. 
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Similarly 


Theorem  2.  Every  modified  Bernstein  function  with  shifted  knots  of  degree  n 
is  a  linear  combination  of  two  modified  Bernstein  functions  with  shifted  knots  of 
degree  n  +  1  : 

G!Uo  w  =  (!4^)  aU«,0  W  +  (^Tf )  oS'w  w  (7) 

where 

<  t  <  and  a,  j3  are  positive  real  numbers  satisfying  0  <  a  <  (3. 

Proof. 

We  have 
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On  using  equalities  (5),  (6),  we  can  easily  get 


^  (^) 


Gkn+1,a,p  (t)+(  —  )Gk+1 


n  +  1 


n-\-l,a,/3 


0 t )•  « 


Theorem  3.  Every  modified  Bernstein  function  with  shifted  knots  of  degree  n 
is  a  linear  combination  of  two  modified  Bernstein  functions  with  shifted  knots  of 
degree  n  —  1  : 

<‘>  -  w.  («) 


where 


<t<  and  a,  ft  are  positive  real  numbers  satisfying  0  <  a  <  (3. 
Proof. 

On  using  Pascal  type  relation  Q'j  =  +  (n^1),  we  get 


(*)  (sis)"  (‘ 


a 


db)  V  »  +  ^/  V»  +  ^ 
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n  —  l\  (n  —  1 

jfe-  1/  +  V  k 
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—  t 
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a  \k[n  +  a  '  n~k 

t - -  - -~t 


V'  n  +  fij  \n  +  /3 


n  —  1  \  1 


+ 


n  —  1\  1 


k 

n  +  (3  (  a 


t  — 


t  — 


a 


n  +  a 
n  +  (3)  \n  +  /3 

k  /  . 

n  +  a 


-  t 


a 


n  +  f3 J  \n  +  (3 


—  x 


n—k 


n—k 


n 


n  +  (3 


)  it)  +  ^  -  tj  W  ( t ). 


When  a  =  f3  =  0,  equalities  (7),  (8)  reduce  to  the  degree  evaluation  formula 
of  the  classical  Bernstein  bases  functions.  ◄ 


3.  Bezier  curves  based  on  modified  Bernstein  bases 

Using  modified  Bernstein  bases  functions  with  shifted  knots,  the  Bezier  curves 
of  degree  n  are  defined  as  follows: 


P(t;a,/?)  =  ^Pk  Gk^(t), 


k= 0 


(9) 
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where  Pk  G  M3  (k  =  0, 1,  •  •  •  ,  n),  and  Pk  are  control  points.  The  adjacent  points 
Pk,  fe  =  0,1,2,  •  •  •  ,n  are  joined  to  obtain  a  polygon  which  is  called  the  control 
polygon  of  Bezier  curves. 

3.1.  Properties 

Theorem  4.  Bezier  curves  based  on  modified  Bernstein  bases  have  the  following 
basic  properties : 

1.  Bezier  curves  have  geometric  and  affine  invariance. 

2.  Bezier  curves  lie  inside  the  convex  hull  of  their  control  polygons. 

3.  The  end-point  interpolation  property:  P(^^;a, /3)  =  Po,  P(^ijrf ;  a, /?)  = 

P„. 

4.  Reducibility:  When  a  =  (3  =  0  formula  (9)  reduces  to  a  classical  Bezier 
curve. 

Proof.  These  properties  of  Bezier  curves  based  on  modified  Bernstein  bases 
can  be  easily  deduced.  ◄ 

Theorem  5.  The  end-point  property  of  derivative-. 
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=  V(t;a,/3 ) 


or 


P  (t;a,/3)  =  V(f;a,/3). 

On  differentiating  both  sides  with  respect  to  ‘t’,  we  have 


Let 


P  '(t;a,P)  =V'(t-a,/3). 


.  ( n\  1  (  a  \k(n  +  a  \  n  k 

Mtia-m=  (*)  (5rr^)b^ 


Then 


V(t-a,P)  =  J2  Pk Ank{t-a,P). 


k= 0 


(A"k(t,a,py=(^j 


n\  1 


'  n  \n 
<  n+{) ' 


k  (  t  — 
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k- 1 


n  +  (3  J  \n  +  (3 


n  +  a 


—  t 


n—k 


n\  1 
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n  \n 
n+/3  / 
n—  1  / 


a  \ 


<-—)  (»-*) 

=  (ra  +  /9){A£lJ(t;a,0)  -  a, /?)}, 


n  +  a 
n  +  /3 


—  t 


n—k—  1 


which  implies 


V'(t;a,P)  =  J2  Pk (A%(P,a,py. 


k= 0 


Now 

V' 

and 

P' 


'  ( p)  =  P'(-^;  a,  /3)  =  (n  +  /3)(PX  -  PoK"  ^  «,  /?) 
\n  +  p  J  n  +  p 


(^:“^)  =  (n+«(p--p»  >(^H 


—  1  +  P  \  n  /  n  —  1  +  a  a 
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Similarly  after  some  computation,  we  have 
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V7 


,  f  n  +  a 
n  +  /3 


;«,/?)  =P  ,(^;«>/3)  =(n  +  /3)(P„-P„-i)A'!:11(^ 


-1  (  n  +  a 


n  +  (5  )' 


W  a,  e)  =  (n+/S)(P„-P„V  A  -  1  +  V»  +  « 


n  +  /3 


ra  —  1 


a 


n  +  /3  n— 1  +  /3 


n— 1 


◄ 


3.2.  Degree  elevation  for  Bezier  curves  with  shifted  knots 

Bezier  curves  with  shifted  knots  have  a  degree  elevation  algorithm  that  is 
similar  to  that  possessed  by  the  classical  Bezier  curves.  Using  the  technique  of 
degree  elevation,  we  can  increase  the  flexibility  of  a  given  curve. 

As  we  know,  Bezier  curves  based  on  modified  Bernstein  bases  with  shifted 
knots  are  given  by 

n 

P(f;a,/3)  =  ^Pk  Gkn^  (t). 

k= 0 

Then  after  applying  degree  elevation  algorithm  on  this,  we  have 

n+ 1 

P(t;a,/3)  =  ]TPk  G£+w  (t) 

k= 0 

as  Bezier  curves  based  on  modified  Bernstein  bases  of  degree  n  +  1,  where 

=  (‘  -  Pk-  +  (sri)  Pk  (12) 

The  equation  above  can  be  derived  from  Theorem  2.  If  we  denote  by  P  = 
(Pq,  Pi,  •  •  •  ,  Pn)T  the  vector  of  control  points  of  the  initial  Bezier  curve  of  degree 
n,  and  by  P(P  =  (Pq  ,  P*,  ■  ■  ■  ,  P*+1)  the  vector  of  control  points  of  the  degree 
elevated  Bezier  curve  of  degree  n  +  1,  then  we  can  represent  the  degree  elevation 
procedure  as 

P(1)  =  Tn+  iP, 


where 
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n  +  1 
n  +  1  —  n 


Tn+ 1  — 


ra  +  1 


0 

0 
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0 
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0 
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0 

0 

n  +  1-2  2 

n  +  1 
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0 

0 


0 

1  1 
n  +  1 


(n+2)x(n+l) 


For  any  l  £  N,  the  vector  of  control  points  of  the  degree  elevated  Bezier  curve 
of  degree  n  +  l  is  P®  =  Tn+i  Tn+ 2  ■  ■  ■  .Tn+ iP.  As  l  — >  00,  the  control  polygon 
P^l  converges  to  a  Bezier  curve. 


3.3.  de  Casteljau  type  algorithm  for  Bezier  curves  with  shifted 
knots 

The  recursive  de  Casteljau  type  algorithm  for  Bezier  curves  with  shifted  knots 
of  degree  n  can  be  expressed  as  follows: 


'  P?(t;a,/3)  =  P?  =  Pi  *  =  0, 1,  2  •  •  •  ,n 

'  pf(«i <*.«  =  (t - pr'(‘;«.« 

r  =  l,  ■■■  ,n,  i  =  0, 1,  2  -  -  -  ,n~  r.,  ^ <t<  0  <a<(3. 

(13) 

Then 

n—  1 

P =  ^P =  •  •  •  =  ^P \{t;a,P)  Gln_rap{t )  =  •  •  •  =  P£  (t;a,/3). 

i=0 

(14) 

It  is  clear  that  this  result  can  be  obtained  from  Theorem  (3).  When  a  =  /3  = 
0,  the  formulas  (13)  and  (14)  recover  the  de  Casteljau  algorithms  for  classical 
Bezier  curves.  Let  P°  =  (Pq,  Pi,--  -  ,  Pn)T  ,  Pr  =  (Pq  ,  P[,  ■  ■  ■  ,  P^_r)r-  Then  de 
Casteljau  algorithm  can  be  expressed  in  matrix  form  as  follows: 

P r(f;  a,  P)  =  Mr(t ;  a,  /?)  •  •  •  Af2(t;  a,  P)Ml(t-  a,  ^jP0,  (15) 

where  Mr(t ;  a,  /3)  is  a  (n  —  r  +  1)  x  [n  —  r  +  2)  matrix  and 


Mr(t;a,P ) 
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0 

0 


n  +  (3 
n 


(m-t) 


0 

0 


(*-3?)  0 

(gg-‘)  (‘-3h  o 

(3?-*)  (*-53)  0 

0  -  (3f -*)  (<-3s)J 


4.  Tensor  product  Bezier  surfaces  with  shifted  knots  on 

I"  a?  n-\-a~ I  y,  I"  ct  n-\-a  ~| 
l_n+/3’  n+/?J  Ln+/3’  n+fi\ 


We  define  a  two-parameter  family  P(w,  v)  of  tensor  product  Bezier  surfaces 
of  degree  m  x  n  as  follows: 

P(rt,  v)  = 


i=0  j=0 


a  n  +  a 

X 

a 

n  +  a 

n  +  (31  n  +  f 

n  +  /?’ 

n  +  (3  _ 

,  (16) 


where  Pjj  G  M3  (i  =  0, 1,  •  •  •  ,m,j  =  0, 1,  •  -  -  ,n),  and  Glm  a  p{u),  G3n  ap(v)  are 
modified  Bernstein  functions,  respectively.  We  refer  to  the  Pjj  as  the  control 
points.  By  joining  up  adjacent  points  in  the  same  row  or  column,  we  obtain  a 
net  which  is  called  the  control  net  of  tensor  product  Bezier  surface. 


4.1.  Properties 

Bezier  surfaces  with  shifted  knots  have  the  following  properties: 

1.  Geometric  invariance  and  affine  invariance  property:  Since 

m  n 

<«>)  =  !>  (17) 

i= 0  j= 0 

P(u,  v)  is  an  affine  combination  of  its  control  points. 

2.  Convex  hull  property:  P (u,v)  is  a  convex  combination  of  P jj  and  lies  in 
the  convex  hull  of  its  control  net. 

3.  Isoparametric  curves  property:  The  isoparametric  curves  v  =  v*  and  u  =  u* 
of  a  tensor  product  Bezier  surface  are  the  Bezier  curves  with  shifted  knots 
of  degree  m  and  degree  n,  respectively.  Namely 


Bezier  Curves  and  Surfaces  based 


15 


u 


e 


a  n  +  a 
n  +  f3'  n  +  (3 


Gm,aj d(V) 


V 


€ 


a  n  +  a 
n  +  j3'  n  +  (3 


The  boundary  curves  of  P  (u,v)  are  evaluated  by  P(u,  P(u, 

p(^,u)andP(^±f,u). 


4.  Corner  point  interpolation  property:  The  corner  control  net  coincides  with 

the  four  corners  of  the  surface.  Namely,  P(^,  =  Po,o,  p(^,  ^+f ) 

p  Tj  (  m+a  \  p  o  /  in+a  n+a  \  p 

r0,n,  rlm+/i>n+a/_r™>  0;  r  f  m+p  ,  n+p  )  ~  rm,n- 

5.  Reducibility :  When  a  =  (3  =  0,  the  formula  (16)  reduces  to  a  classical 
tensor  product  Bezier  patch. 


4.2.  Degree  elevation  and  de  Casteljau  algorithm 

Let  P  (u,v)  be  a  tensor  product  Bezier  surface  with  shifted  knots  of  degree 
m  x  n.  As  an  example,  let  us  consider  obtaining  the  same  surface  as  a  surface  of 
degree  (m+  1)  x  (n  + 1).  Hence,  we  need  to  find  new  control  points  P*  ■  such  that 


m  n  m+1  n+1 

p(u,„)  =  =  EE  P  i,jGm+  l,a,  ^U)Gn+l,a,^V)- 

i= 0  j= 0  i= 0  j= 0 

(18) 

Letai  =  l-^,  /3j  =  l-n±l=l 

Then 


Ptj  =  ai 


Pj  Pi-i,j-i  +  oti  (1  —  /3j)  Pi-i, j  +  (1  —  ai)  (1  —  f3j)  P ij, 


(19) 


which  can  be  written  in  a  matrix  form  as 


[m+l— i] 
[m+1] 


[m+l— i] 
[m+1] 


Pj— 1,7  —  1 

.  PiJ-1 


1  - 


[n+l-j]' 


n+1] 
[n+l-j] 
[n+1] 


De  Casteljau  algorithm  can  also  be  easily  extended  to  evaluate  points  on 
a  Bezier  surface.  Given  the  control  points  Pt,j  6  M3,i  =  0,1,---  ,m,  j  = 
0, 1,  •  •  •  ,w,  we  have 
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P v)  =  P P  j0  -  Pijt  =  0,  X,  2  •  •  •  ,  m;  j  =  0, 1,  2  •  •  ■ 


P- /(«,«)  = 


n+P 


m+p 


-t 


n+P 


t~ 


m+p 


jr  — l,r  — 1 
*+ 1,3 


r  =  1, 


,  k,  k  =  min(m,  n)  i  =  0, 1, 2  •  •  •  ,  m  —  r;  j  =  0, 1 , 


■  i,j+ 1 

^f- l,r—  1 

■  i+l,j  +  l  . 

,  n  —  r, 


ra+/3  /  n+a  _  j.\ 
n  Vn+/3  *7 

(f  _  _2_\ 

n  \L  n+p) 

'  Q:  <  +  <C  n+a 
n+p  -  1  —  n+p  * 
(20) 


When  m  =  n,  one  can  directly  use  the  algorithms  above  to  get  a  point  on 
the  surface.  When  m  /  n,  to  get  a  point  on  the  surface  after  k  applications  of 
formula  (20),  we  use  formula  (15)  for  obtaining  the  intermediate  point 

Note:  We  get  classical  Bezier  surfaces  for  (u,v)  G  [^qrg,  ^|]  x 
when  we  set  the  parameter  a  =  f5  =  0. 


5.  Concluding  remarks 

Bezier  curves  and  surfaces  are  constructed  with  the  help  of  modified  Bernstein 
bases  functions  with  shifted  knots  for  t  G  [^qrg>  ^+f]  •  Introduction  of  parameters 
a  and  (3  enable  us  to  shift  Bernstein  bases  functions  over  subintervals  of  [0, 1]. 
We  get  Bernstein  functions  on  [0, 1],  when  we  set  the  parameters  a  =  f 3  =  0. 

Simulation  study  is  performed  through  MATLAB  R2010a.  It  has  been  found 
out  that  Bezier  curves  that  are  generated  over  any  subinterval  of  [0, 1]  based  on 
modified  Bernstein  bases  functions  which  satisfy  the  properties  of  non  negativity, 
partition  of  unity  and  end  point  interpolation  property  are  similar  to  the  Bezier 
curves  that  are  generated  based  on  classical  Bernstein  bases  functions  over  the 
interval  [0,1].  This  result  is  shown  in  Figure  3  where  Bezier  curves  based  on 
modified  Bernstein  bases  functions  is  shown  by  blue  curve  get  overwrite  by  Bezier 
curves  based  on  classical  Bernstein  bases  functions  shown  by  green  curve. 


6.  MATLAB  Code 


m= [4,0]  ; 

1= [5 , 0]  ; 

cpx=  [37  40  39  29  23  26  45  ] 
cpy=  [38  37  27  26  36  50  56  ] 


°/0  value  of  alpha 
l  value  of  beta 

l  x-coordinate  of  control  points 
l  y-coordinate  of  control  points 
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(a)  (b) 


(c) 

Figure  3:  Bezier  curves  based  on  modified  Bernstein  bases 


for  shifloop=l:2 

i=i; 

a  =  m(shifloop);  f3  =  l(shifloop); 

n=6 ; 
d=l; 


for  y  =  :  0.01  :  ?±§ 


n+/3 


n+/3 
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d=d+l ; 
end 


bezierx= [1 :  d-1] ; 
beziery= [1 : d-1] ; 


bezierbernsteinx=0 ; bezierbernsteiny=0 ; 


for  k=0:n 


end 

bezierx{i)  =  bezierbernsteinx  * 
beziery(i)  =  bezierbernsteiny  * 

i=i+l ; 
end 

if  (shif loop==l) 

c=plot (bezierx , beziery) 
set(c, 'Color’ , ’blue’ , ’LineWidth’ ,2) 

else 

c=plot (bezierx, beziery, ’g’) 

set(c, ’Color’ , ’green’ , ’LineWidth’ ,2) 
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end 

hold  on 

end 

plot ( cpx , cpy , ’ — k* ’ ) ; 

hlegl  =  legend ( ’Modified  Bernstein  bases’ , ’Classical  Bernstein  bases  (green 
overwrite  blue) ’,’ control  polygon’ ,’ location’ , ’northwest ’) ; 
set (hlegl , ’Box’ , ’off ’ ) 
hold  on 
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